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We study single-qubit gates performed via stimulated Raman adiabatic passage (STIRAP) on 
a spin qubit implemented in a quantum dot system in the presence of phonons. We analyze the 
interplay of various kinds of errors resulting from the carrier-phonon interaction (including also 
the effects of spin-orbit coupling) as well as from quantum jumps related to nonadiabaticity and 
calculate the fidelity as a function of the pulse parameters. We give quantitative estimates for an 
InAs/GaAs system and identify the parameter values for which the error is considerably minimized, 
even to values below 10 -4 per operation. 
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I. INTRODUCTION 

Quantum dots (QDs), among many other systems^ 
are considered to be promising candidates for an imple- 
mentation of quantum information processing schemes. 
Due to their atomic-like structure^ one can easily sin- 
gle out a subset of states to encode the logical qubit 
values. In principle, these systems provide for sta- 
ble coherent memory if the information is encoded into 
the long-living electron spin^ which motivated a spin- 
based proposal for quantum information storage and 
processing^ On the other hand, experimental demonstra- 
tions of coherent control over the charge (orbital) degrees 
of freedom 5 - 6 - 7 8 - 9 and the recently performed two-qubit 
gate based on a confined biexciton systemifi prove the 
feasibility of quantum coherent manipulation of carrier 
states on picosecond time scales. It has been therefore 
proposed^ to implement the qubit states as the vacuum 
and single exciton states in a QD, switched by resonant 
optical coupling and providing the two-cubit conditional 
gating via inter-QD dipole-dipole interaction. 

Both the spin-based and the charge-based proposals 
suffer from serious difficulties. The spin switching time 
in typical structures is very long due to weak magnetic 
coupling. The orbital degrees of freedom do not provide 
for long operation times due to the finite exciton lifetime, 
usually of order of 1 nsiiSii^ It seems therefore natural to 
seek for a scheme in which the logical values are stored 
using spin states, while the operations are performed via 
optical coupling to the charge degrees of freedomi^i^, 
also using QD systems in QED cavitiesiiSii 7 . 

A promising solution, proposed recently^ is to encode 
the qubit states into spin states of an excess electron in 
a QD and perform operations by employing the stim- 
ulated Raman adiabatic passage (STIRAP) to a state 
localized spatially in a neighboring dotii^ (An alterna- 
tive scheme not relying on the auxiliary state has also 
been proposed^). The STIRAP technique uses three 
laser fields that couple the two qubit states as well as 
the auxiliary state to a fourth state, a charged exciton 
(X~ , or trion), composed of two electrons with oppo- 



site spins and a hole. In the presence of laser fields with 
slowly varying amplitudes, the system evolves adiabat- 
ically, following the states of the interacting system of 
carriers and electromagnetic field (trapped states). The 
driving fields may be chosen in such a way that the trion 
state is never occupied (in the ideal case) so that the 
scheme is not affected by the decoherence resulting from 
its finite lifetime. It can be shownSi that with prop- 
erly chosen phases of the laser pulses a pre-defined qubit 
superposition gets coupled and undergoes an adiabatic 
transition to the second dot and back which shifts its 
phase by a desired angle with respect to the other, or- 
thogonal superposition that remains decoupled from the 
laser fields. This results in an arbitrary rotation of the 
qubit state around an arbitrary axis on the qubit Bloch 
sphere. 

The essential difference between the atomic systems, 
where such quantum-optical schemes are successfully 
applied^ and the solid state QD systems, where their 
new implementation is proposed, is the nature of the 
environment. In high-quality samples at low tempera- 
tures the dominant coupling to the external degrees of 
freedom is that involving lattice modes (phonons). The 
coupling mechanisms include interaction with lattice po- 
larization (longitudinal optical, LO, phonons) and with 
piezoelectric fields induced by phonon-related strain (lon- 
gitudinal and transverse acoustic, LA and TA, phonons) 
as well as the effective influence of strain-induced band 
shift, described in terms of the deformation potential 
coupling to LA phonons. Even restricted to acoustic 
phonons, this kind of external bath shows various pe- 
culiarities compared to models usually assumed in gen- 
eral studies^ Its low-frequency behavior depends on the 
coupling mechanism and on the wave-function geome- 
try and is always super-ohmic, i.e., its spectral density 
grows super-linearly with frequency^ Due to the local- 
ization of carrier wave-functions on a scale much larger 
than the lattice constant, a high-frequency exponential 
cut-off in the effective phonon spectral densities appears 
well below the Debye frequency. Moreover, apart from 
the non-diagonal coupling terms describing real transi- 
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tions, there is usually a diagonal coupling which leads 
to pure dephasing effects 2 ^* 2 ^ resulting from the lattice 
relaxation after a fast (compared to phonon frequencies) 
change of the carrier statei 2 ^* 2 ^ Such an effect manifests 
itself in optical experiments as a fast partial decay of the 
signal coherenc e!?' 12 in excellent agreement with theoret- 
ical modeling assuming its phonon-related origin. 27 - 28 

The characteristic time scales of these intrinsically non- 
Markovian pure dephasing processes are determined by 
the localization (QD size) and are typically much shorter 
than any real phonon- induced transition process. More 
importantly, they overlap with the time scales proposed 
for optical qubit control^ It has been shownSi that the 
demand to avoid these pure dephasing effects limits from 
below the gating times, thus shrinking the time scale win- 
dow defined, on the other side, by the long-time decohcr- 
ence processes (e.g. the exciton lifetime). 

In this paper we study the influence of the coupling to 
the phonon degrees of freedom on the fidelity of the single 
qubit rotation via the STIRAP process 2 -^ implemented in 
the double-QD structured Even if the possible phonon- 
assisted transitions to other states may be neglected, the 
diagonal terms still give rise not only to pure dephas- 
ing effects but also to transitions between the trapped 
carrier-field states. The probability of these phonon- 
induced transitions becomes very high if the spacing be- 
tween the trapped energy levels falls into the area of 
high phonon spectral density and the overall error is 
roughly proportional to the process duration. Such high 
error rates are critical for quantum computation schemes 
where extremely high fidelity is required (e.g., errors not 
higher than ~ 10~ 4 per gate are allowed for two-qubit 
operations) in order to provide for scalable devices in- 
cluding quantum error correction schemes. We discuss 
how these strong decoherence processes may be avoided 
by either decreasing the trapped level separation (low- 
frequency regime, exploiting the super-ohmic behavior 
of spectral densities) or increasing it beyond the cut-off 
(high frequency regime). We show that in both cases 
one encounters a trade-off situation, due to the oppo- 
site requirements for phonon-induced jumps (short dura- 
tion) and for the fundamental adiabaticity condition and 
pure dephasing (slow operation): In the low- frequency 
regime, avoiding phonon-induced transitions contradicts 
the condition for avoiding nonadiabatic jumps between 
the trapped states, which may be overcome only by con- 
siderably extending the process duration. In the high- 
frequency case, there is a competition between the pure 
dephasing and the phonon-induced transitions that is 
overcome by increasing the trapped state splitting, tak- 
ing advantage of the particular structure of the phonon 
spectral density for a double dot structure. 

The paper is organized as follows. In the next Sec- 
tion [H] we present the general derivation of the phonon- 
induced error for an arbitrary system evolution. The 
Section lnTl describes the model of the specific system dis- 
cussed in the paper and derives the carrier-phonon cou- 
pling relevant for our discussion. The Section lTvl provides 



a description of the STIRAP qubit rotation procedure for 
completeness and necessary reference. In the central Sec- 
tion |V| the results of Sec. El are applied to the STIRAP 
procedure described in Sec. II VI with the phonon perturba- 
tion derived in Sec. IIIII This section contains also some 
general discussion. In the Section IVH we present the re- 
sults for specific pulse shapes in order to get some quan- 
titative estimates for an InAs/GaAs QD system. Finally, 
the Section IVIII summarizes and concludes the paper. In 
addition, some technical details and further analysis, in- 
cluding the effect of the spin-orbit coupling, are presented 
in the Appendices. 

II. PHONON-INDUCED DECOHERENCE: 
GENERAL THEORY 

Subject of this paper is the optically induced dynamics 
in a quantum dot structure coupled to a phonon bath. In 
this Section we derive the equations for the reduced den- 
sity matrix of the carrier subsystem in the leading order 
in the phonon coupling, assuming that the unperturbed 
(ideal) evolution of the noninteracting system, described 
by the unitary evolution operator 

U {t)=U c {t)®e- iH ^ h (1) 

is known (see also Ref. l30|) . Here, Uq is the evolution 
operator for the carrier subsystem coupled to the external 
light field in absence of carrier-phonon interaction and 
-ffph is the free phonon Hamiltonian. The relevant carrier 
states in the quantum dot are assumed to form a discrete 
set |n), n = 0, 1, 2, . . ., and the phonons are described by 
destruction and creation operators bk and 6^ referring to 
bulk phonon modes with wave vector k. 

The interaction between the carriers and the phonon 
modes is written in the general form 

V = J2 S nn>®Rnn>, (2) 

nti' 

where S nn i = S n , n = \n)(n'\ act in the Hilbert space of 
the carrier subsystem and 

Run' = R\, n = ]T F nn , (k) (b k + &!„) , (3) 
k 

with F nn >(k) = F* t (—k), affect only the phonon envi- 
ronment. 

We assume that at the initial time — to the system is 
in the product state 

£>(-*o) = Po ® Pr, Po = l^oXV'ol, (4) 

where \i/jq) is a certain state of the carrier subsystem 
and pr is the thermal equilibrium distribution of phonon 
modes. Physically, this is justified by the existence of two 
distinct time scales: the long one for the carrier decoher- 
ence (e.g. about 1 ns ground state exciton lifetim e) 12 ! 13 ) 
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and the short one for the reservoir relaxation (1 ps pure 
dephasing timeA2*2ii2ii) . 

The starting point is the evolution equation for the 
density matrix of the total system in the interaction pic- 
ture with respect to the externally driven evolution {To, 
in the second order Born approximation with respect to 
the carrier-phonon interaction 31 



Q(t) = e(-t ) + I dr[V(T),Q(-t )] (5) 



/ dT dr'[V(r), [V(r% Q (-t )]}, 
n J-t J-to 

where 

g(t) = Ut(t) Q (t)U (t), V(t) = ul(t)VU Q {t). 
The reduced density matrix of the carrier subsystem is 

p(t) = U c (t)p(t)Ul(t), p{t) = [Tr R ~g(t)} , 

where the trace is taken over the reservoir degrees of free- 
dom. Note that in this paper the symbol g always refers 
to a density matrix in the full carrier-phonon Hilbcrt 
space while p refers to reduced density matrices either 
in the phonon or the carrier subspace. The first (zeroth 
order) term in J^J obviously yields 

P M(t) = u c (Wo)(Mul(t) = \Mt))(Mt)l (6) 

The second term vanishes, since it contains the thermal 
average of an odd number of phonon operators. The 
third (second order) term describes the leading phonon 
correction to the dynamics of the carrier subsystem, 



(J) 



dr / dT'Tr R [V(T),[V(r'),g(-t )]]. 



pW(t) 



The first of the four terms resulting from expanding 
the commutators in J7J) is —QtPa, where 



1 



0* - s^z. 



nn' mm 



dr / dr' 

to J — to 



X S nn ' (r)S mm i (t )(Rnn'{ T T )Rmm')- 

The operators S and R are transformed into the interac- 
tion picture in the usual way 

S nn '(t) = U%(t)S m >U (t), R nn ,{t) = U^(t)R nn ,U Q (t) 

and (O) = Ti r [Opt) denotes the thermal average (obvi- 
ously [Ua(t), pt] — 0). Using the symmetry of the oper- 
ators S nn i and R nn > the second term may be written as 



— poQt- I n a similar manner, the two other terms may 
be combined to <!>t [po]: where 



nn mm 



xS>n n '(T )pS mrn > (t) (Rrnm' (j~ T )Rnn') ■ 

In terms of the new hermitian operators 

A t = Q t + Ql, ht = ^(Q t - Ql), (8) 

the perturbation to the density matrix at the final time 
t J7J may be written as 

p (2 Ht) = -i [ht,/*] - \ {A uPo } + * t \po]. (9) 

The first term is a Hamiltonian correction which does not 
lead to irreversible effects and in principle may be com- 
pensated for by an appropriate modification of the con- 
trol Hamiltonian Hq . The other two terms describe pro- 
cesses of entangling the system with the reservoir, leading 
to the loss of coherence of the carrier state. 

Introducing the spectral density of the reservoir, 



Rnn' : mm' (^) 



one may write 



1 



2wh 2 



dt{R nn ,{t)R mm ,)e^\ (10) 



nn' ,mm' 



nn mm 



(ii) 

where the frequency-dependent operators have been in- 
troduced, 



Y nn i(u}) 



S rm ,(T)e luJT dT. 



(12) 



Using 1)10(1 one has also 

Q t = I duj [ dT [ dT'Q{T-T') 

i i J J —tn J — tn 

nn' mm' u 11 

x Snn' {j~)Smm' (t )Rnn' ,mm' ( w ) e ^ • 

Next, representing the Heaviside function as 



B(t) = -e luJt J 



2iri to' - u + i0+ ' 



we write 
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Qt = 



/ dujR nn ^ rmn >{ 

nn' mm' 

5^ / didR nn ' 7 mm'(> 
nn' mm' 



dLo' Y^, n {J)Y mml {uj') 
2iri u)' - oj + iO+ 

rh / 



-m8{J - lu) +V — 



where V denotes the principal value. 

Hence, the two Hermitian operators defined in (JSJ take 
the form 

A t = ^2 E J dcjR nn \ mm '(uj)Y^ n (uj)Y mm/ (uj) (13) 



nn mm 



and 



EE 

nn' mm' 



dcjR 



nn'mm 



(14) 

du/y^ n (u/)F w (u/) 



2n 



In the following, we will be interested in the system state 
at the final time t = +t.Q, after all the pulses have been 
switched off. 

In the quantum information processing context it is 
customary to quantify the quality of the operation in 
terms of the fidelity, which is a measure of the overlap 
between the desired (unperturbed) state and the actual 
final state, F = Tr[U c (t)paUl(t)p(t)}. The error is then 
defined as the fidelity loss, 5 = 1 — F. From Eqs. JJJJ and 
ijffil one has 

5 = -<Vo|/5 (2) |^o> (15) 

^] J d^J Rnri ,mm' (w) 
nn'mm' 

x(V>oK t ,^)P ± Y mm ,(o;)|V>o>, 

where P 1 - is the projector on the orthogonal complement 
of |V>o) in the carrier space. In this order the unitary 
correction generated by ht does not contribute to the 
error. 

The calculation presented above requires two input 
components: the specific form of the interaction poten- 
tial [Eq. J5J] for a given problem and the unperturbed 
time evolution [Eq. QJ]. These two necessary elements 
are derived for our qubit system in the two following sec- 
tions. 



III. THE QUBIT SYSTEM AND ITS 
INTERACTION WITH PHONONS 

In the following part of the paper, the general theory 
will be applied to a specific system of two quantum dots 
containing one excess electron and coupled to the trion 



state in order to perform an arbitrary rotation in the 
qubit space by means of the STIRAP. Here we formu- 
late the model of this system and derive the Hamiltonian 
describing its interaction with the phonon environment. 

The Hamiltonian describing this system and its cou- 
pling to lattice modes may be written as 



H = H C + flp h + V, 



(16) 



The first term is the STIRAP Hamiltonian including 
both the qubit states and the control fields. The 
implementation 18 defines the qubit by two a y spin eigen- 
states of a single excess electron in one of the QDs 
("large") from a vertically stacked pair. In order to 
perform a general single-qubit rotation between the two 
qubit states |0) and |1) an auxiliary state |2) is used2i, 
in which the electron is shifted to the second ("small") 
dot and has the same spin orientation as in |0). All these 
three states are coupled to a fourth state |3), a charged 
exciton (trion) state, by laser beams fio,^i,^2, respec- 
tively. The Hamiltonian for such a system in rotating 
wave approximation (RWA) is 

h c = E e >X"i 

n 

2 

+ ^,nn n (t) (e^-'-^lnXSl+H.c.) , 

n=0 

where e„ are the energies of the corresponding states, 
the slowly varying pulse envelopes O n (£) are real and 
positive, u> n are the corresponding laser frequencies and 
and S n are the phases of the pulses. This Hamiltonian 
induces the unitary evolution described in the previous 
section by the operator Uq [Eq. |QJ]. 

In order to achieve the Raman coupling, the frequen- 
cies uj n of the laser beams must be chosen such that the 
detunings from the corresponding dipole transition en- 
ergies £3 — e„ are the same for all the three couplings. 
Therefore, we put u> n = e^/h — e n /% — A, n = 0, 1, 2, 
where A, the common detuning, is one of the parameters 
to be tuned for optimal performance. In the rotating 
frame, defined by \h) = e^'-^ln), n = 0,1,2, the 
RWA Hamiltonian i|17|) may be written 



(18) 



7iA|3)(3| + nn n (t)(e-^\h){3\ + e^\3)(h\), 



n=0 
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where S n — 5 n — So (only the relative phase of the pulses 
matters). 

The second term describes the free phonon evolution, 



H, 



where /3l, (3k are phonon creation and annihilation opera- 
tors (with respect to the crystal ground state). Through- 
out the paper, the phonon branch index will be implicit 
in fc, unless it is explicitly written. Together with Hq 
[Eq. (|17|l ]. the above phonon Hamiltonian describes the 
known, unperturbed evolution of the system, given by 
Eq. <|TJ>- 

The final term is the carrier-phonon interaction. Since 
the adiabaticity inherent in the STIRAP procedure ex- 
cludes the possibility of inducing high-frequency dynam- 
ics and also all the trapped state splittings should be 
at most of several meV (to avoid crossing with excited 
carrier states) , the discussion will be restricted to acous- 
tic phonons. The Hamiltonian describing the electron- 
phonon interaction in the coordinate representation is 



(19) 



where r denotes the electron coordinate (a similar con- 
tribution appears for holes). The coupling constants 
for the longitudinal and transverse phonon branches are, 
respectivel y 32 ^ 3 , 



,0) - 



(20a) 



2p c V n LJi(k) 



de 

ak - i Mi(fe) 



and 



v 



(ti,t 2 ) 



(20b) 



de 



2p c V n u t (k) £ £ s 



M HM (k)F nn ,(k), 



where i,ti,t 2 refer to the longitudinal and two transverse 
acoustic phonon branches. Here e denotes the electron 
charge, p c is the crystal density, V n is the normaliza- 
tion volume of the phonon system, 0Ji,t are the phonon 
frequencies, d is the piezoelectric constant, £o is the vac- 
uum dielectric constant, £ s is the static relative dielectric 
constant and a is the deformation potential constant for 
electrons. The functions M s depend on the orientation of 
the phonon wave vector—. For the zinc-blende structure 
they are given by 



M s {k) = 

2 k x ky{^& s k) z -\- kyk z [&s,k^)x H~ k z k x ^s,k)y i 



(21) 



where k — k/k and e Si fc is the unit polarization vector 
for the wave vector k and polarization s. 



In the basis of the confined states relevant for the STI- 
RAP process the carrier-phonon interaction Hamiltonian 
(Unj reads 



V= \n)(n'\Y,fnrAk)((3 k + (3l 



n.n'—O 



k ; j 



(22) 



where, for single-electron states, f nn '(k) = VkJ- n n'(k) 
with the form factors T nn i(k) depending on the wave 
function geometry and given by 

^nn'(fc) = / d*r** n (r)e ikr * n ,(r), (23) 

where ^ n (r) is the envelope wave- function of the elec- 
tron. The coupling constants f nn >(k) include all the cou- 
pling mechanisms relevant for a given phonon branch and 
have the symmetry f nn '(k) = /*<„(-&)■ 

We will assume that the two spin states used to encode 
|0) and |1) correspond to the same orbital wave- functions 
so that the couplings foo(k) and /n(fe) are equal. The 
couplings /oi(fc), /io(fc), /i2(fc) and / 2 i(fc) vanish since 
the spin orientation in the state |1) differs from that in 
|0) and 1 2) (the effects of the spin-orbit coupling are dis- 
cussed separately below). Moreover, it is assumed that 
there is no overlap of wave-functions between the states 
|0) and 1 2), so that also f 02(h) and f 20(h) vanish. 

An important point is that, since the electron resides 
normally in the large dot, at the initial moment the lat- 
tice is relaxed to the corresponding minimum ( "dressing" 
of the electron in the coherent deformation field). This 
may be accounted for by defining the modes with respect 
to this shifted equilibrium, so that the ground state of 
the interacting system corresponds to the new phonon 
vacuum, i.e., by transforming to new phonon operators 
bk according to 



/oo(fc) 



(24) 



Upon transformation to these new modes the interaction 
reads 



V = E \n)(n\J2Fnn(k){b k + bl k ) 

n=2,3 k 
' 2 

Ei n X 3 iE/™ 3 ( fc ) ( b ^+ b -k) +h.. 



.n=0 



where F nn (k) = f nn (k) — foo(h)- Moreover, the carrier 
Hamiltonian undergoes a renormalization which is, how- 
ever, inessential for our discussion. In the rotating frame 
the above interaction Hamiltonian reads 



n=2,3 k 
' 2 

5>)<3|5> n3 (fc) (b k + bl k )+R, 



.71=0 
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where F n3 (k) = / n 3(fc) e ~ i(wrat ~' 5o) - This Hamiltonian 
is of the form (J2J with S nn i = and R nn > = 

J2k^nn'(k)(b k + 6_ fe ). The spectral densities [Eq. (fTTHl ] 
have the explicit form 

Rrm',mm'(u) = ^2 XI ( k ) F ™>m ( k ) ( 2 6) 

x + l)(5(w - Wfe) + nfe<5(w + Wfe)] , 

where are phonon occupation numbers. Note that 
not all possible couplings appearing in the general form 
of Eq. J2J are present in our case. It is clear from Eq. 
(|25|l that the phonons influence the dynamics only when a 
transfer from the large dot to a spatially different carrier 
state (small dot or trion state) occurs. 

The interaction potential given by Eq. 1(25(1 will be used 
in the calculation of phonon-induced decoherence accord- 
ing to the general theory of Sec. [H] First, however, one 
has to describe the unperturbed evolution which is the 
second necessary ingredient of the calculation. This is 
done in the following section. 



Thus, the parameters x an d Si define two orthogonal 
states in the qubit space. The laser pulses affect only one 
of these states, the coupled (bright) state \B), while the 
other orthogonal combination, \D), remains unaffected. 

At a fixed time t, the Hamiltonian 127fl has the eigen- 
statcs 

I do) = cos(9|B) — e-^ 2 sin6»|2) (28a) 

|o_) = cos</>(e^ 2 8inO\B)+ cos 6\2)) (28b) 

-e lS2 sin0|3), 

\a+) = sm(/){sm9\B}+e- tS2 cos6\2)) (28c) 
+ cos0|3), 

where 

tan o = — — , sm o = —= 1 , 



IV. THE STIRAP PROCEDURE FOR A 
SINGLE-QUBIT ROTATION 

In this section we present the formal descriptions^*^ of 
the stimulated Raman adiabatic passage without exter- 
nal perturbation. Along with the results of the previous 
section this will allow us to use the general theory of Sec. 
[n]for the description of phonon-induced dephasing. 

The system is modeled by the Hamiltonian given by 
Eq. 118(1 . The envelopes of the first two pulses, Oo^, are 
chosen proportional to each other so that they may be 
written as 

Sl (t) = tt Q1 (t) cosx, fli(t) = fl 01 (t) sinx, 

with a certain parameter x S (0; §) defining the fixed 
ratio of the pulse intensities. In terms of the new basis 
states 

\B) = cos X |0)+e-^sin X |l), 
\D) = -sinx|6) +e~ J,5l cosx|l), 
the Hamiltonian (|18|) now reads 

H C = ftA|3><3| + ~iM*)(|£><3| + |3><S|) (27) 

+ ~fi 2 (i) (e-^|2)(3|+e^|3)<2|). 
I 



The corresponding eigenvalues are SAo,±, where 

x = o, A± .WE+H±M, (29) 

The system evolution is realized by an adiabatic change 
of the pulse amplitudes (see Fig. in this application, 
the detuning remains constant). Initially (at the time 
—to), both pulses are switched off, hence 4> = 0, then fl 2 
is switched on first, hence also 6 = 0. Therefore, |ao) co- 
incides with \B) and |o_) with |2). During an adiabatic 
evolution of the parameters, the states move along the 
corresponding spectral branches. During the first pas- 
sage, 82 = and 9 is changed from to tt/2. At the end 
of this stage, when the pulses are switched off (<fr = 0), 
the electron is in the state — 12). The second passage 
takes 8 back from tt/2 to 0. Now, however, 62 7^ so 
that the adiabatically followed system state is e tS2 \ao) 
and the final state is e lS2 \B). Note that the desired sys- 
tem evolution relies on the angle 9 determined by the 
ratio ( so cane d mixing angle), while the absolute 

value of these pulse amplitudes remains a free parameter 
that may be used for optimization against decoherence 
effects. 

The evolution operator corresponding to this proce- 
dure may be written (in the basis \B), |2),|3)) 



Uc(t) 



e lS2 cos ( 



sm ( 




e -tA_+t& cos d> sin ( 



_ e -t\_+ts 2 gin( 



e + sm <p sm 6 



(30) 
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FIG. 1: (a) An example of pulse shapes (solid) and the result- 
ing structure of the dresses levels (dashed) . The arrows show 
the phonon-assisted transitions, as described in Sec. [3 (1) 
the pure dephasing effect, (2,3) the transitions between the 
trapped states, (b) The evolution of the functions 9 (dashed) 
and sin 28 (solid) for the pulse sequence shown in (a) . 



where 9,4>, X± are slowly varying functions of time, 62 = 
for t < and 



A±(t) 



/ drA±(r). 
Jo 



As shown in Ref. l2ll the phase shift of the bright state 
resulting from the procedure described above is equiva- 
lent to the rotation in the qubit space |0), |1) around the 
axis determined by \ an d by the relative phase <5i be- 
tween fio and Qi. The rotation angle is equal to the 
62 phase of the second pulse sequence. The character- 
istic feature of the STIRAP is that no special form of 
the pulses is required. Thus, from the point of view of 
the unperturbed evolution, the detuning A and the pulse 
envelopes f2oi,2 are to a large extent arbitrary. This free- 
dom may be used for minimizing the perturbing effects 
of the environment. 

Ideally, the state |2) is only occupied during gating, 
while the state |3) is never occupied. This is true un- 
der the assumption that the evolution is perfectly adia- 
batic. However, any change of parameters can never be 
infinitely slow and the probability of a jump from |ao) to 
one of the two other states \a±) remains finite. In the 
lowest order, the corresponding probability amplitudes 



,34 



c±(t) 



dT(a±(r)\ip(T)) exp 



lj\ ± (r')dT' 



where ip(t) is the state evolving adiabatically from the 
initial one. Let us write the general initial state in the 



^o)=co S ||B)- e ^sin||D). (31) 

The qubit rotational is performed by two well sepa- 
rated, mirror-symmetric pulse sequences differing only 
by a phase. Thus, using the explicit formulas JSSaJc) one 
may write 

c±(t) = cos | e -*A±(*) [c±-c* ± ], (32) 



where 



c± = 



dr 



sin 0(r) 
cos (J)(t) 



6(7 



JA±(t) 



(33) 



with sin</>(r), cos0(t) corresponding to '+' and '— ', re- 
spectively, and the integral involves only one pulse se- 
quence. If the evolution induced by the pulse sequence 
is symmetric with respect to a certain time t\ (the time 
around which the pulse sequence arrives) the amplitude 
(|33|l may be written in the form c± — ie* A± ' tl ^|c±|. 

In order to discuss the general properties of the nona- 
diabatic jump amplitudes, let us write the evolution of 9 
in the form 



9(t) = 



TO 



(34) 



where 9 is a function of unit width, so that To is the time 
scale of the evolution of 9. (The total duration of the 
gate, including two pulse sequences, is roughly an order 
of magnitude longer). If the functions <fi(t) and A±(i) 
change slowly around t = t\, then 



sin 4>{ti) 
cas(f>(t\) 



0o[roA±(ti)], 



where 



du9'( 



is a function of unit width, with a fixed value at x = and 
vanishing for x 3> 1 (here prime denotes the derivative 
with respect to the argument u). Hence, the nonadiabatic 
jump amplitudes l|32|) are small when 



|A±| > l/r , 



(35) 



which is the standard adiabaticity condition. 

It is interesting to note that for symmetric pulses the 
final transition probabilities (|32() vanish for A±(ti) = 
(n + 1/2)tt. This fact is due to destructive interference 
of the jump amplitudes during the first and the second 
pulse sequence. Although it might be tempting to ex- 
ploit this cancellation and perform a successful passage 
for times and laser beam parameters that do not satisfy 
the condition (|35|l . such a procedure requires a detailed 
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knowledge of the excitonic dipole moments and a pre- 
cise control over the laser beam properties. Moreover, 
the cancellation takes place only in the final state, while 
during the process the other states are occupied, which 
leads to the non- vanishing occupation of the X~ state 
and to decoherence, contrary to the original motivation 
of this qubit implementation. In order to avoid these ef- 
fects, the envelopes of the transition probabilities should 
be used as the actual bound to the nonadiabatic-jump- 
related error. 



V. INTERACTION WITH THE PHONON BATH 
DURING THE STIRAP PROCESS IN A QD 
SYSTEM 

In this section we apply the general theory from Sec- 
tion [n] to the qubit rotation performed via a STIRAP 
process, as described in Section llVI implemented in the 
double- QD system. 

In the Hamiltonian (|25|l . the only non- vanishing non- 
diagonal coupling is F n3 (k). Let us note, however, that 
for this coupling one has, according to (|12[) . 



)t+s„ 



Uc„ m Uc3 m '\m)(m'\, 



For the initial state (|31|l . using the explicit evolution 
operator (|30l) . the contributions from the three states 
\4> n ) =srn|i?|B)+e i ^cos|'!?|£»),|2) ) |3) are, respectively, 

p CO 

Sl (ui) = - - sin d / dte~ wt sin 2 6{t) , (38) 



«2,3(w) 



W 1? 
" 7T COS — 

2 2 



cos 4>{t) 
sm.<f)(t) 



die 



sin20(i)f 



(39) 



-iA=F (t) — 162 



These three contributions correspond to transitions indi- 
cated graphically in Fig. ^ 

Following the argument that led to Eq. 1)32(1 . these 
functions may be written in the form 

\s n (uj)\ = 2u n (0)Re[si(w)], 

where u\ = isin?9, 1*2,3 = cos(i?/2), and 



dte 



'■ sin20(i)0(i) 



(40a) 



ie-^lSiHI, 



where Ucnm are the elements of the evolution operator 
(|30|l , varying at most with frequencies ~ X± . It is there- 
fore clear that this function is peaked around uj « -w B , 
i.e., at the optical frequencies which are many orders of 
magnitude higher than any phonon frequencies present in 
Rnn',mm'(w) [Eq. I|10fl]. Thus, inter-band non-diagonal 
phonon couplings do not contribute to l|15[) . This is con- 
sistent with the rotating wave approximation and may 
also be understood by noting that the second Born ap- 
proximation accounts for processes that may be repre- 
sented as a series of emission and absorption processes 
involving arbitrarily many photons but only one phonon. 
Each photon process takes the system from the states 
0,1,2 to 3 with the exchange of a large energy while 
a non-diagonal phonon process produces the same state 
change but with negligible energy exchange. Thus, en- 
ergy can never be conserved in a process involving the 
inter-band phonon term. 

Since the adiabatic evolution Uc does not transfer 
qubit states into |3), U^\3) remains orthogonal to \B). 
Hence, 133(0;) does not contribute to (|15|l and we may 
write 



dio 



R(u) 



S(uj), R(u) = i?22,22(w), 



(36) 



with 



5h = w 2 ^k^i4miv^)i 2 = Ei s «hi 2 ' ( 3? ) 

n n 

where the sum runs over a complete set of states \tp n ) 
orthogonal to \iffo). 



« 2|3 (w) = -\\ dte-^ t+A ^-^ 



COS (j)(t) 

sin cj)(t) 



-i[u>t 1 +A T (t 1 )+8 2 /2]i~ 



sin 20(i) 
(40b) 



where the integrals are now over one pulse sequence and 
the final equalities hold for symmetric pulse sequences. 

Using the representation l|34() of the system evolution 
and denoting the Fourier transform of 8' sin 20 by gi (x) 
we find |si(a>)| = <7i(wto). Since t\ ^> tq and |si(w)| 2 is 
integrated with the slowly varying spectral density, the 
oscillating terms do not contribute and one may write 

| Sl H| 2 « isin^lsx^To)! 2 . 

Hence, the function, s\(uj) is centered at lu — and 
broadened by a factor 1/tq due to the time-dependence. 
It is responsible for the pure dephasing effect m& The re- 
sulting error, according to (|37)|) . will grow with broad- 
ening of si(w), i.e., with decreasing process duration. 
Hence, similarly to the fundamental condition l|35|l . it 
always favors slow operation. However, it is independent 
of the trapped level splittings and reflects only the low- 
frequency properties of the spectral density (at a given 
temperature). For R(ui) w Ro^ n , n > 3 this pure de- 
phasing error at the temperature T is 



j(pd) 



-(n-1) 



k B T ~(n-2) 
T 1 



ksT » %/tq. 



(41) 



It should be noted that the crossover from the low to 
high temperature behavior is governed only by the pulse 
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duration (irrespective of the system parameters) and for 
durations of the order of 10 ps it takes place at T ~ 0.1 
K. 

By a similar argument, the two other functions may 
be approximately written as 



l*2,3M| S 



(wro) 2 



i) 



cos 
sin 2 



g\ [(w + A t )t ] , 



where g<i{p) is the Fourier transform of sin20(u). These 
functions have a similar 1/tq broadening but are also 
shifted to the spectral position u> = — A±. They de- 
scribe the error resulting from phonon-assisted transi- 
tions between the trapped states |flo,±) (see Appendix 
for further support to this interpretation). In view 
of the condition 135|) . this shift must be larger than the 
broadening and for rough estimates the latter may be ne- 
glected (if the spectral density varies slowly on the scale 
of this broadening; the role of oscillations in the spec- 
tral density is discussed below). Hence, one may write 

6 (tr) = j(tr 



6^'\ where 



, R(u) , l2 



(42) 



R(-X T ) J du 



32,3 M 



(J 



The error is therefore proportional to 
R(\±)t . 



(43) 



Thus, for a fixed spectrum of the trapped states this error 
grows linearly in time (in the leading order), which is a 
usual characteristics of real transition processes. 

In order to maximize the fidelity of the coherent op- 
eration, one must find the trade-off between the errors 
caused by phonon-assisted transitions which favors 
short process durations, and the other two restrictions, 
related to pure dephasing (|41(l and to the general adi- 
abaticity condition (|35f) both increasing for fast evolu- 
tion. As can be seen from the orders of magnitude of the 
spectral characteristics determining the error [Fig. [21 and 
Eq. J3BJ], in general, the fidelity may be strongly de- 
creased. However, contrary to the simple excitonic qubit 
case^S, the STIRAP procedure in a QD system provides 
two ways to avoid this limitations. 

First, due to the super-ohmic properties of the phonon 
spectral density, R(ui) ~ w™, n > 2, all error contri- 
butions may in principle be minimized by locating the 
trapped levels X± in the low-frequency sector and de- 
creasing them while simultaneously increasing the gate 
duration To. 

Second, the values of A± may be chosen sufficiently far 
beyond the cut-off frequency. The contribution from the 
phonon-induced transitions and nonadiabaticity effects 
may then be arbitrarily small and the error is limited by 
the pure dephasing effect, restricting the possible gate 
speed-up. However, one should keep in mind that in the 




2x10 



1x10 



FIG. 2: (a) The functions |s n (u;)| 2 describing the phonon- 
induced errors (for pulses as in Fig. and the total spectral 
density of the phonon reservoir R(ui) at T = (solid) and 
T = 5 K (dashed) for the model InAs/GaAs system (Tab. |TJ. 
The inset shows the exact shape of one of the spectral features, 
(b) The contributions to the spectral density at T=0: DP 
coupling to LA phonons (solid) and piezoelectric coupling to 
TA (dashed) and LA (dash-dotted) phonons. Inset: high- 
frequency behavior with the two bounds defined in Appendix 

m 



high frequency domain there may be additional reservoir 
excitations (including two-phonon processes) that are not 
accounted for in this model. 

The error effects discussed here originate from the in- 
teraction between phonons and the orbital degrees of 
freedom used to operate the qubit. On the other hand, 
for a spin qubit one expects some contribution to the de- 
coherence induced by the spin-orbit (SO) coupling. The 
electron confined in a quantum dot does not interact with 
other carriers so that, in contrast to highcr-dimcnsional 
systematise, dephasing of the electron spin requires 
interaction with phonons^ or nuclei^. In the Appendix 
IP we analyze the former channel. We show there that 
the Markovian decay of spin states in our system is 
very slow and leads to negligible error over the times 
relevant for the qubit operations. On the other hand, 
non-Markovian SO-related effects induce the transitions 
between the same states as the direct phonon coupling 
but are many order of magnitude weaker due to very 
small SO-induccd phonon coupling resulting from rela- 
tively large energy separation of orbital states. Thus, the 
SO-related effects do not affect the discussion presented 
here. 
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Electron effective mass 


* 

m 


0.067m c 


Static dielectric constant 


e B 


13.2 


Piezoelectric constant 


d 


0.16 C/m 2 


Longitudinal sound speed 


Cl 


5600 m/s 


Transverse sound speed 


Ct 


2800 m/s 


Deformation potential for electrons 


(7 


-8.0 eV 


Density 


Pc 


5360 kg/m 3 


Lande factor 


9 


-0.44 


Spin-orbit coupling constants: 






Rashba 


a 





Dresselhaus 


P 


1 nm/ps 


Level separation 




46 meV 


Electron wave-function widths: 






in-plane 


l± 


5.0 nm 


z-direction 


lz 


1.5 nm 


Dot separation 


D 


6.0 nm 



TABLE I: The GaAs material parameters an d QD system 
parameters used in the calculations (after Refs. I40ll4 11) . 



VI. QUANTITATIVE RESULTS FOR A MODEL 
PULSE SEQUENCE 

In this section we calculate the errors for a STIRAP 
operation on a single qubit performed with specific pulse 
shapes. In order to get quantitative estimates and to 
identify the key error-inducing mechanisms in various 
regimes of operation we use the material parameters and 
QD characteristics for an InAs/GaAs system which is 
frequently used as the "typical" system for the proposed 
qubit implementations. The system parameters are col- 
lected in the Tabled 

It is known that the STIRAP procedure is rather in- 
sensitive to the exact pulse shape. In order to simplify 
the discussion, we choose the pulse sequence 



-[(t±t 1 )/r ]2 



1/2 



which results in a very simple form for the time depen- 
dence of the mixing angle, 



sin 26 = 



1 



'*±*i 



The envelope 51 onv (t) may be any function approximately 
constant around t\. For the numerical calculations we 
take 



^env(^) 



1 



1 + a cosh 



with a = 10~ 4 , t\ — 0.4r (Fig. ^ corresponds to this 
pulse choice), fl 2 is proportional to the total power of 
the three pulses. The constant fi, along with A, must be 
tuned for minimizing the decoherence effect. 



For such a pulse sequence one finds explicitly from Eqs. 

up 



ft 

|c±| 2 = 47TCOS 2 - 



< 47rcos 2 



sin 
cos 2 



sin 



sin 2 [A±(ii)]e 



-2(A±t ) 2 



-2(A±t ) 



where the envelope of the oscillations has been taken as 
the safe bound to the error, in accordance with the dis- 
cussion in Sec. IIVI For the purpose of analytical esti- 
mates the values of </> = <j>{ti) and A± = \±(t%) are as- 
sumed constant. The resulting error is equal to the sum 
of the two transition probabilities |c±| and depends on 
the initial state l|31|l . since c± = c±(i?, ip). The error 
averaged over the initial states is 



(5( na ) — 
1 



(44) 



— I dtp I ^sintf(|c + (tf,^| 2 + | C _(^)| 2 ) 



2tt sin 2 <^r 2(A - ro) " + cos 2 0e" 2(A+ro) 



The spectral functions Si(uj) relevant for the phonon- 
induced dephasing are 



\ Si (lj)\ 2 = ^sin 2 t?— — sin^ile-*^ 2 (45) 

2 2 + 7T 

and 



|s 2:3 H| 2 = 2 7 rcos 2 -(o;r ) 2 cos 2 M 1 +A T (t 1 )](46) 



cos" 



-(u;+A t ) 2 t 2 



The total error is calculated as the sum of the nona- 
diabatic jump probability l|44f) and the phonon-induced 
contributions given by (|3b137fl with the spectral functions 
(I45|l and l|46|). The phonon spectral density correspond- 
ing to our model double-dot InAs/GaAs system is derived 
and discussed in the Appendix lAl and plotted in Fig. [21 

The resulting error, averaged over (?9, ip) as in Eq. 144(1 . 
as a function of the pulse intensity parameter and de- 
tuning A for a fixed process duration tq is shown in Fig. 
The nontrivial interplay of the three error contribu- 
tions discussed above together with the oscillating high- 
frequency tail of the phonon density of states R(ui) (see 
inset in Fig. |5p) lead to an intricate parameter depen- 
dence of the total error. There are clearly several pa- 
rameter combinations for which the error becomes small. 
With the help of the formulas l|29|) one finds that the area 
(0a) corresponds to A_ in the low-frequency region, while 
in (0b) A+ is small and A- shifted beyond the phonon cut- 
off. The valleys (1), (2) . . . correspond to A_ positioned 
at one of the minima in the high-frequency tail of R(lu) 
and A+ shifted beyond the thermal cut-off for phonon- 
absorption processes, i.e. hX+ > k^T. For T = the 
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■10 -5 

A[ps-i] 



FIG. 3: The dependence of log 10 8 on the pulse parameters A 
and f2 at T = and T = 1 K for to = 50 ps. Numbers refer 
to the parameter regimes discussed in the text. 



absorption processes are not allowed at all and these ar- 
eas are not separated from the (Ob) region. 

The detailed analysis of the error value along the (0a) 
valley at various temperatures (Fig. shows that at 
T ^ the dependence is not monotonous. The absolute 
minimum always corresponds to very low fl and A, for 
which both trapped states X± lie in the low frequency re- 
gion. At high frequencies, the error values reach a plateau 
after passing (at T > 0) through a second, very shallow 
minimum (due to the subtle interplay of the error contri- 
butions weighted by the parameter-dependent sine/) and 
cos(/> factors). In between, there is either a monotonous 
increase (at T — > 0) or a transition through a local max- 
imum, as the A+ state crosses the frequency sector with 
high spectral density for phonon absorption (cf. Fig|5J|. 
Fig. shows the interplay between different error con- 
tributions when the Rabi frequency f2 is changed for a 
fixed detuning A. In this range of parameters, for the 
specific system under study, the pure dephasing contribu- 
tion turns out to be small compared to the errors related 
to real phonon-induced transitions and to nonadiabatic 
jumps which create a trade-off situation with one or two 
well-defined parameter sets corresponding to the minimal 
errors. 

The above results show that for a fixed pulse duration 
tq the error values are bounded from below, precluding 
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t = 50 ps 
A= 1.2 ps -1 
T= 1 K 
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FIG. 4: (a) The dependence of the error for growing detuning 
with A_ = const, along the (0a) minimum in Fig. [3] (b) The 
total error (solid) and the individual contributions: nonadia- 
batic jumps (dashed), phonon-assisted real transitions (dot- 
ted) and pure dephasing (dash-dotted) for a section of the 
parameter space. 



a perfect operation for any parameter values. However, 
due to the super-ohmic behavior of all the contributions 
to the phonon spectral density (at low frequencies), the 
total error is decreased when the process time grows and 
the trapped level splittings decrease. The minimum er- 
ror achievable for different process durations at various 
temperatures is plotted in Fig. [5Ji,b and the correspond- 
ing laser beam parameters are shown in Fig. [SJ:,d. Both 
the values at the global minimum (Fig. |5p,,c) and at the 
shallow local minimum (Fig. \Sjp, d) are shown. In or- 
der to allow for any subtle interplay of parameters, for 
each To the full minimization with respect to both A and 
O was performed. As expected, the error decreases for 
longer pulse durations, but the decrease is only polyno- 
mial (S ~ 1/to at higher temperatures and To > 10 ps). 
Therefore, rather long pulse durations are necessary to 
reduce the error considerably. Moreover, the optimiza- 
tion is obtained for rather unusually small parameter val- 
ues (Fig. \5]p) and is very sensitive to their precision. Still 
another restriction is that in this low-frequency regime 
the optimum is searched against the nonadiabatic jump 
error and is reached for ToA± > 1. As soon as To becomes 
comparable to the trion radiative lifetime (~ 1 ns), the 
optimal value of A± falls within the broadening of the |3) 
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FIG. 5: (a,b) The minimal achievable error as a function of 
the process duration corresponding to the optimal pulse pa- 
rameters with both A± in the low-frequency area (a) and with 
A+ in the high frequency area (b) (for T = the plateau value 
for A — + oo is shown). (c,d) The optimal pulse parameters 
(detuning and Rabi frequency) realizing the minimal error for 
these two configurations. The legend in (d) applies to both 
(c,d). 



state, disabling the adiabatic passage. 

The parameter dependence of the error in the (0b) area 
is in a way analogous. Here, however, it is |A_| that 
must be shifted far beyond the positive frequency cut-off. 
Even at zero temperature, the positive-frequency part of 
the spectral density extends to relatively high frequencies 
(with oscillations manifesting themselves as local minima 
in Fig. Oh)- Therefore, this parameter regime is always 
less favorable than the previous one. 

In view of the limited possibility of fidelity optimiza- 
tion in the low-frequency region for reasonable process 
durations, it is interesting to study the high-frequency 
parameter range. In contrast to the previous case, the 
values of fl and A may now seem unusually high, but 
the results of Fig. [3] show that by increasing the split- 
ting between the trapped state energies the error may in 
principle be reduced to arbitrarily low values. 

The Figure HH shows the error along the (1) and (2) 
areas (Fig. |3J) for fixed pulse duration at various temper- 
atures, as well as the contributions to the error in one 
case. The trapped states are now split by several meV, 
so that the nonadiabatic error is negligible (except for 
sub-picosecond pulses). However, the speeding-up of the 
dynamics is limited by the pure-dephasing contribution. 
On the other hand, extending the pulse duration is un- 
favorable due to the phonon-assisted transitions. The 
interplay of these two contributions for a given pulse du- 
ration, temperature, and A_ yields a series of minima, 
corresponding to A + traveling across the oscillations of 
R(u>), as shown in Fig. [fj] Note that at low temperatures 
only one minimum exists, belonging actually to the (0b) 
parameter area, but at higher temperatures the absolute 




Q [meV] CI [meV] 

FIG. 6: The error as a function of Q for A_ = const, along the 
1 (a) and 2 (b) areas in Fig. [3] In (a) the individual contri- 
butions to the error at T — 5 K are also shown: nonadiabatic 
jumps (dashed) , phonon-assisted real transitions (dotted) and 
pure dephasing (dash-dotted). 

minimum shifts to the high-frequency region. 

The minimum value reached depends on the pulse du- 
ration, with a certain optimal trade-off which depends, 
however, on the chosen value of A_ and decreases sub- 
stantially for subsequent minima of the spectral density. 
The resulting minimum value, obtained by numerical 
minimization with respect to fl and A for a range of pulse 
durations, is shown in Fig. |JJi,b. The individual contri- 
butions shown in Fig. EH show that pure dephasing in- 
deed limits the fidelity for short pulses but in the optimal 
duration range the non-monotonous To-dependence of the 
error is determined exclusively by the phonon-assisted 
transition contribution. This astonishing effect is in fact 
due to the relatively narrow minimum of R(u>) in which 
|A_| is placed. For short pulses, S3 (a;) becomes broad 
(pure dephasing broadening of the A_ level), increasing 
the overlap with R(u>). For large r , the linear increase 
of S due to long process duration becomes dominating, 
leading to a minimum at a certain point. 

VII. CONCLUSIONS 

We have studied the fidelity of the coherent operation 
on a QD spin qubit rotated by a stimulated Raman adi- 
abatic passage to a neighboring dot and back. We have 
shown that, in addition to the usual limitation of the 
speed of an adiabatic process, the presence of the phonon 
reservoir imposes two further restrictions: The transfer 
must be slow in order to minimize the pure dephasing 
effect but it should not take too long in order to avoid 
transitions between the trapped carrier-light states. The 
general formalism was applied to an InAs/GaAs self- 
assembled system of typical size. It turns out that for 
most values of pulse parameters (pulse intensities and de- 
tuning) in meV range the error is high enough to totally 
prevent the coherent operation. However, there are also 
narrow parameter areas where the fidelity is considerably 
higher. 

The super-ohmic characteristics of the spectral density 
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FIG. 7: (a,b) The minimal achievable error as a function of 
the process duration corresponding to the optimal pulse pa- 
rameters with both \± in the high-frequency parameter areas 
1 (a) and 2 (b) and the optimal pulse parameters (detun- 
ing and Rabi frequency) realizing the minimal error for these 
two configurations (c,d). In (a) the contributions to the er- 
ror at T = 1 K are shown: nonadiabatic jumps (dashed), 
phonon-assisted real transitions (dotted) and pure dephasing 
(dash-dotted). 



associated with the phonon reservoir admits minimiza- 
tion of the total error by increasing the duration of the 
process while simultaneously decreasing the trapped level 
energies. However, the pulse durations necessary for a 
considerable reduction of the error in this low-frequency 
regime arc of order of hundreds of picoseconds which 
leads to nanosecond overall gate durations (full sequence 
of two pairs of pulses). Moreover, the resulting trapped 
state energies become extremely small, approaching the 
typical lifetime broadening of the trion state used for the 
Raman coupling. 

It is found that the qubit operation may be per- 
formed with much higher fidelity if the trapped states 
are pushed beyond the cut-off of the effectively coupled 
phonon modes. An additional advantage comes from the 
oscillatory structure of the phonon spectral density for 
a double-dot system. In this way the error at T = 
may be reduced to the value of ~ 10 -3 and well below 
10~ 4 for the trapped state energy splitting of 4 meV and 
8 meV, respectively (for the system geometry assumed 
here). The latter values lie in the spectral region where 
the acoustic phonon effects dominate the decoherence, 
well below any spectral features (LO phonons, higher ex- 
citon states) not included in the discussion. It is remark- 
able that such low error values are achieved with pulse 
durations of the order of 10 ps which, compared to the 
long electron spin decoherence time, even up to tens of 
milisecondssS,, opens a broad time window for a large 
number of gating operations. 

We have analyzed also the errors related to the spin- 
orbit coupling. It turns out that these effects are negli- 



gible in the present implementation. 

These optimistic conclusions are somewhat limited by 
the strong temperature dependence of phonon occupa- 
tions, especially in the low-frequency regime, leading to 
a fast increase of the error at non-zero temperatures. In- 
deed, in some cases the minimal error may grow even by 
an order of magnitude as soon as the temperature reaches 
1 K. 

The strong dependence of the phonon-related error 
on the material parameters and system geometry opens 
some possibility of system engineering and optimiza- 
tion. For example, the high-frequency asymptotics of the 
phonon spectral density is governed by the QD height: 
higher dots assure a faster decay. On the other hand, 
in the low-frequency sector the phonon spectral density 
scales with the square of the inter-dot distance, favoring 
rather flat structures. Also increasing the lateral size re- 
duces the phonon coupling but, at the same time, lowers 
the excited states restricting the high-frequency range of 
operation. This shows that finding the optimum may 
be nontrivial and may depend on the frequency sector 
chosen for the qubit operation. It should be noted that 
the high-frequency spectral density is dominated by the 
deformation potential coupling which is present in any 
semiconductor system but in the low-frequency domain 
the piezoelectric effects dominate. This might suggest 
using non-piezoelectric materials. 

Let us note also that the single-qubit error calculated 
in this paper gives also an estimate of the two-qubit op- 
eration if the latter is performed using dipole coupling 
between the auxiliary states in the STIRAP scheme. 18 
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APPENDIX A: PHONON COUPLINGS AND 
SPECTRAL DENSITIES 

In this Appendix we derive the spectral density of the 
phonon reservoir i?22(w) and study its properties for low 
and high frequencies. 

The phonon coupling constants F-22(k) = f 22(h) — 
/oo(fe) have the same structure as the original constants 
[Eq. I|20al b)] with the form factor replaced by F(k) — 
F12 (&) — Foo(k). Let ^-L,s(fe) denote the form factors, 
calculated according to Eq. I|23|) . for the ground-state 
electronic wave-function in the large (L) and small (S) 
dot. Assuming that the dots are stacked along the z axis 
at the distance D, one has 

T{k) = e l£ ^T s {k) - e- iS ^F L (k). 
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The long- wavelength properties of the coupling constants 
do not depend on the wave-function geometry. Indeed, 
Fs,h(k) = 1 + 0(k 2 ) and T{k) = iDk z + 0(fc 3 ). 

The coupling constants for arbitrary k depend obvi- 
ously on the specific form of the wave-functions. For 
simplicity, we assume Gaussian wave-functions, 



^ , L,s( r ') = Nexp 
Then 



2 \ I 



r± 



_L.S 



'zL,S 



(Al) 



Allowing for a small difference between the dot sizes we 
write Zi LjS =ll± ±A(?i), l 2 zLS =l 2 z ± ±A(l 2 z ), so that 



2i sin 



-m 2 



Dk z , klA(ll) + k 2 A(ll) 



Dk z 



cos ■ 



Hence, the size difference brings only a small correction 
and will be neglected. 

Assuming isotropic phonon dispersions, the spectral 
density R(ui) = i?22,22(w) [Eq. Q1U|)] may be written as 

RM = ^E/« 2 

x [(n k + l)S(u) - w fc ) + n k 6{uj + w fc )] 
x-^ J cos9d8 J dip\F$(k)\ 2 . 

The LA phonons are coupled both via piezoelectric and 
deformation potential interaction. However, due to dif- 
ferent inversion symmetry the mixed terms vanish upon 
angle integration and the two terms contribute indepen- 
dently. 

The deformation potential term is 



i?( DP >M = 4 DI V[nBH + i]/ (DP) M, 



where 



R 



(DP) 



1 D 2 oi 



3(2tt) 2 hp c cj 
and the function /( DP ) (u>) is defined as 

r ir/2 



/ (DP) M = £ 



d8 cos ( 



tt/2 



4 sin 2 ( sin I 



(A2) 



1 /u)l± 
xexp -- 

2 \ ci 

so that J( DP ' (w) -> 1 as w -> 



cos 



I 2 

1 ± 



sm 



For the piezoelectric contributions we choose the 
phonon polarizations 

ei.fc = fe = (cos 6 cos 4>, cos 6 sin sin 6) , 
e t i.fe = (- sin0, cos0, 0), 
et2.fc = (sin cos <j>, sin sin ^, — cos 9); 

then the functions M s [Eq. I)2ip] are 
3 

Mi = — sin 20 cos sin 2</>, 

Mti = sin 29 cos 20, 

Mt 2 = (3 sin 2 6< - 1) cos sin 2(j>. 

The corresponding terms in the spectral density are 
R^(u) = Y^<- S) ^[n B {u) + l]/( p -)(c), (A3) 



where 



R 



1 /edDY 
2fip c (27r)3cf /ifl J 



i /-Tr/2 4 sin I ^p- sm f 
/ (p ' s) H = — / d0cos0M a 2 (0) 

Ms J-tt/2 ( Du^ ~ 



3i n 2 (^sine) 



2ci 



2 i 2 

cos 0+7#- sin ( 



/( p > s )( w ) -> 1 as w 0, and 

M 2 {6) = I d<pM 2 (9,cp), n s = I decos9sm 2 6M 2 (6). 



The specific values are /xj = ^ = 167r/35, /xt 2 = 
167r/105. Thus, the low-frequency behavior of the in- 
dividual contributions to the spectral density is ~ w 3 
and r~j lo 5 for the piezoelectric and deformation potential 
coupling, respectively. 

The behavior in the high frequency limit is deter- 
mined by the coupling to phonons with wave vectors 
in the strongest confinement direction, i.e. along the 
z axis. The piezoelectric coupling in this direction is 
suppressed by the geometrical factors M s and the corre- 
sponding contributions to the spectral function decrease 
rapidly. Moreover, the frequencies of TA phonons are rel- 
atively low and the piezoelectric coupling to LA phonons 
is much weaker. The frequencies of LA phonons reach 
much higher values, e.g. over 20 meV for GaAs, and 
their dispersion remains approximately isotropic and lin- 
ear up to several meV^ 1 Expanding the integral into an 
asymptotic series one finds an upper estimate for l|A2fl . 



/( DP) 



< 



12c? 



D 2 (l 2 ± -l 2 )cu 



1 e -*M 
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In vicinity of the points u> n = 4mrc\/D, the following 
lower bound approximately holds 



/ (DP) H 



> 



3cf 



1 e zi'ci) 



(A4) 



(see Fig. |2t>) The oscillatory behavior of the spectral 
density for large frequencies follows from the fact that 
the predominant contribution in this sector comes from 
phonons along the strongest confinement direction, lead- 
ing to a pronounced destructive interference of interac- 
tion amplitudes in the double-dot structure aligned along 
this direction. 



APPENDIX B: TRANSITIONS BETWEEN THE 
TRAPPED STATES: FERMI GOLDEN RULE 

In this Appendix we show that the error 5^ [Eq. (|42|l ] 
may be interpreted, in terms of the Fermi Golden Rule 
(FGR) , as resulting from transitions between the trapped 
states of the confined electron in the external driving field 
[transitions (2) and (3) in Fig.^. 

Inserting the definition H40b|) into the Eq. (|4^|l and 
performing the frequency integration (for X± , <f> w const) 
we get 



cos 2 ^i?(-A T ) 



J dt sin 2 26{t). (Bl) 



Let us now consider the probability of phonon absorp- 
tion or emission leading to a transition from the state 
\ao) to \a±). The duration of a single absorption or 
emission process is of the order of the inverse phonon 
frequency (i.e. trapped level spacing). Hence, in view of 
the adiabaticity condition l|35|) this process is fast com- 
pared to the characteristic time scale of the system evo- 
lution. Therefore, it is reasonable to calculate the FGR 
probability for absorption or emission at fixed values of 
system parameters and include the time-dependence re- 
lated to the STIRAP passage only at the level of the 
rate equations. Assuming the initial state (|3ip. taking 
the matrix element of the phonon coupling Hamiltonian 
(|25|l between the trapped states J35a|-c) and applying the 
FGR in the standard form one finds for the transition 
probability 



2tt 1 



i) 



. - cos — 

h 4 2 



cos" 

. 2 

sin 



sin 2 2^)]T|F(fc)| ; 



[d(h\ T - huj k )n k + 6(h\ T + hu> k )(n k + 1)] 



7T , 

2 2 



sm 



sin 2 28(t)R(-X T ). 



Solving the rate equation for the jump probability with 
the above time-dependent rate w(t) we find the error 
probability for the whole process duration 



^(ti) _ i _ ex p 



-_(t)dt 



(B2) 



For small error values this reduces to (|B1|) . However, it 
gives also an estimate for the error beyond the applica- 
bility of the perturbative treatment. 



APPENDIX C: SPIN-FLIP EFFECTS DUE TO 
SPIN-ORBIT COUPLING 

In the present Appendix, we discuss the additional er- 
ror due to the presence of spin-orbit coupling for the elec- 
tron. We will show that each of the spin-conserving de- 
phasing channels discussed in the main body of the paper 
is accompanied by a spin-flip channel which is, however, 
several orders of magnitude weaker in a self-assembled 
system. There is, moreover, an additional error related 
to a spin-flip transition in the small dot but it is also 
extremely small. 

We start the quantitative analysis by adding the spin- 
orbit coupling to the qubit Hamiltonian in Eq. <|16[) . 



Hr 



H 



so 



H L (t). 



Here H d = p 2 /(2m*) + U(r), where m* is the electron 
effective mass and U(r) is the confinement potential; 
Hz = (l/2)g^B<7 y is the Zeeman energy (g is the ef- 
fective Lande factor, fiB is the Bohr magneton and a y is 
the Pauli matrix; the magnetic field is oriented along y), 
Hi,(t) describes the coupling to the control laser field and 
Hso = f3(~p x Vx+PyVy)+a(PxVy-PyVx) is the spin-orbit 
term composed of the Rashba and Dresselhaus coupling 
with the constants a and f3, respectively. 

Following Ref. HE we l°°k f° r the unitary transforma- 
tion e s that eliminates the spin-orbit coupling from the 
stationary Hamiltonian Hi = H^ + Hz + -ffso- To the 
leading order in the SO coupling one has 

e s H ie ~ s = H d + Hz + H so + [S, H d + Hz]. 

For the harmonic confinement U{r) = (1 / '2)m* <jJq(x 2 + 
y 2 ) + (l/2)m*w 2 z 2 , uj z » lu , the SO coupling is pertur- 
batively eliminated with the choice 



S 



. 9PbB 

1 (M 2 



(f3p x - ap y ) 



where we omitted an irrelevant position-dependent part. 

Upon the canonical transformation, the electron- 
phonon Hamiltonian l|19(l becomes, in the leading order 
in the SO coupling, V = V + V a , where the additional 
term is 



V„ - [S,V] 
~ \huj) 2 



where v k are defined by Eqs. f2()al b). Within the reduced 
subspace spanned by the relevant states, this operator 
has non-vanishing elements only between those states 
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that have overlapping wave functions and opposite cr y 
spins. Thus, one has 



[F^ S) (k)(\0)(l\e- lEzt ^ +H.c.) 

fc 

+^ L )(fc)(|2)(2'| e -^ zt /' l + H.c.)l {fih + Plu), 



where Ez = gHB,B is the Zeeman energy splitting, |2') is 
the state in the "small" dot with flipped spin and 



F^ L \k) = F( s < L >(-k 



(CI) 



h(/3k x - ak y )v k Ts^{k)e 



±- * 
t- 



with the form factors given by Eq. (lAlT) . 

Upon the phonon equilibrium shift given by Eq. Q24[l. 
the above interaction Hamiltonian produces a small spin- 
dependent renormalization of the qubit Hamiltonian. 
More importantly, the canonical transformation implic- 
itly performs a transition to the eigenstates of the full 
Hamiltonian including the SO term. These states may 
couple to the control field in a different manner than the 
original states, which is reflected in the present formalism 
by the correction terms resulting from the transformation 
H L (t) = e s H L (t)e- s = H L (t) + (t) + . . .. The fea- 
sibility of the ST1RAP process in the presence of such 
spin-dependent terms is a separate problem, far beyond 
the scope of the present paper. Here we assume that the 
control field can be appropriately modified so that the 
new states may be evolved according to the same STI- 
RAP transfer as the original ones. The weakness of the 
spin-dependent phonon effects, as discussed below, sug- 
gests that also these SO corrections might be of minor 
importance. 

The SO contribution to the error may be written, in 
analogy to Eqs. (|15I36[) . as 



5 a = I (k}R a (w)S a {u), 



where 



with 



(C2) 



S„(w) =Lu 2 J2\(iPn\YAu)\4> }\ 2 , (C3) 



Y a (cu) = / dte lut ul(t) 

J ~t 

x(|0)(l|e- ijEzt/?i + \2)(2'\e- lEzt/h + R.c.)U c (t). 



The summation in Eq. (|C3|I now involves all the states 
\ip n ) orthogonal to \ipo), including |2'). 

Similarly as in the Appendix ^ we find the low- 
frequency expressions for the spectral densities [Eq. (fTUfl ] 
corresponding to the SO coupling (|C1|) via the two dif- 
ferent coupling channels 

4 DP >(a,) = i^% 5 [n B M + l], (C4) 



where 



and 



where 



R 



(dp) _ l {gmBy ha 2 2 2 

-° - 12^ (^ )4 p c cj [a +Ph 



4 p ''^) = J&M»bM + i], (C5) 



tH (hujo) A p c c° s \e e s J 

with 71 = 1/35, jt x = 7t 2 = 2/105. In the frequency 
range typical for the Zeeman energies in GaAs at mod- 
erate magnetic fields (~ 0.1 meV) the piezoelectric cou- 
pling to transverse modes dominates. 

In the high-frequency region, where the deformation 
potential coupling dominates, we find the asymptotic es- 
timate 



(C6) 



4tt 2 (Hluo) 4 p c cf 
xKMll]e-» ( t» i ^ , 

The spectral function S a (u>), pertaining to the driven 
evolution of the system, depends in a complicated way on 
the performed qubit rotation and on the initial state. In 
order to reduce the complexity, we restrict the discussion 
to a 7r/2 qubit rotation around the a x axis, i.e. x = 7r /4, 
<5i = 0. We parameterize the general initial qubit state in 
the form (|31|l and calculate the spin-flip contributions to 
the error averaged over the Bloch sphere of initial states. 
We will restrict the discussion to the most interesting 
high-frequency regime of operation, where the following 
hierarchy of time scales may be assumed: to ^> fi/Ez 3> 
h > r . 

Under these simplifying assumptions and neglecting 
terms proportional to EzTq the first contribution to the 
spectral function S a (uj) [Eq. (|C3|l ] is 

l<o|^M|i)| 2 « 

— {2t - 2t x ) [5{lo + Ez/H) + S(uj - Ez/%)] 



2 

av ' 



where 'av' denotes averaging over the Bloch sphere of 
initial states, as in Eq. I|44|). The first contribution leads 
to the Markovian (Fermi Golden Rule) spin-flip proba- 
bility over the time 2to — 2t\ during which the electron 
is located in the first ("large") dot, up to the factor 1/3 
resulting from averaging of the spin-flip rates for various 
superposition states. The probability of such a process 
in a self-assembled quantum dot is extremely low due 
to large confinement energy. Indeed, the spin-flip rate 
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pertaining to this contribution is, according to Eq. 1C2I1 . 1 [cf. Fig. [3], i.e. A_ = 12 ps 1 , one finds 



2tt 

w = —[B„{E z /h) + R a {-E z /ri)} 



1 16 ( wbB \ 
3 105tt V Kojo J 

1.0 ■ 10" 



h p c c 2 



de 



(C7) 

(« 2 + /? 2 ) 



where we used the low-frequency formula (|C5I) for the 
spectral density (the piezoelectric coupling to transverse 
phonons dominates in this sector), substituted the Zee- 
man energy E z — g[i B B and assumed that k B T > E z . 
The final value corresponds to T = 1 K and B = 1 T. 

The other term describes an additional contribution to 
the spin-flip transition closely related to the pure dephas- 
ing effect [Eq. (|41|) ] . Since the low- frequency behavior of 
the spectral densities for direct and SO-induced processes 
[Eqs. (US! and l|U5f) ] is the same, the ratio between the 
SO effect and the direct pure dephasing is 



5P_ 

<5(pd) 



R 



(P,t) 

<r0 



h 2 (a 2 +[3 2 )( gf i B B) 2 



R, 



(P,t) 



D 2 



1.7-10" 



where we included again only the dominating contri- 
bution from the piezoelectric coupling to transverse 
phonons. It is clear that the SO-related process is negli- 
gible compared to the pure dephasing. 

Under the same assumptions as above, we have for the 
two other SO-induced contributions 

K2,3|y(w)hto>l2v = Als2, 3 HI 2 . 

As discussed in Sec. the functions S2,3(w) are rela- 
tively sharply peaked around to = — A T . Since the latter 
values are large, we use the high-frequency asymptotics 
for the spectral density (now the deformation potential 
coupling to longitudinal phonons dominates) and the ra- 
tio between the SO-related process and the transitions 
described in Sec. Ivlmav be written as 



^ 3) 

5(2,3) 



id DP) (-A T ) 
i?( DP )(-A T ) 

(W?) 2 h 2 



D 2 c 2 T 



\i«i-H){* 2 +n 



where we used the asymptotic formula (|A4(I since the 
optimal pulse parameters correspond to Azp in a local 
minimum of the spectral density. For the parameter area 



5]? 



1.8-10 



-in 



Again, the SO-related process is many orders of magni- 
tude weaker than the transition discussed in the paper. 

Apart from these contributions there is another one, to 
the spin- flipped state |2') in the small dot. This process 
is not possible in the absence of the SO coupling. The 
relevant spectral function is 



l<2'|l»l^>l 



sin 2 [(u- E z /h)ti+g(u)} 



(w - E z /ny 



h 2 {w-E z /h), 



where we denoted 



dt cos 66e l 



ig(uj) 



h(uj). 



One has for the STIRAP transfer g(0) = 0, h(Q) = -1. 
The width of h(uj) is of the order of l/r . For large t\ 
(long dwelling time in the small dot), this yields the spin- 
flip transition rate according to the Fermi Golden Rule, 
analogous to Eq. (|C7fl up to an averaging-related factor. 
For typical dwelling times ~ 100 ps this would produce 
a negligible contribution to the error of order of 10~ 14 
for B ~ 1 T, T ~ 1 K. However, under the assumptions 
made above the FGR is not applicable; instead one may 
approximate 

|<2'|yM|Vo)l 2 v « \h 2 {u). 

Using the low-frequency and high-temperature approxi- 
mation to R^'^^lo) one finds 



K a0 



B<jO 



k B T 1 

h T 
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-14 



at B ~ 1 T, T ~ 1 K. Note that under the conditions 
assumed here, the dominating contribution comes from 
the dynamical effect characterized by the inverse pulse 
duration 1/ro. 

In conclusion, the Markovian spin-flip rate for a self- 
assembled QD with typical level separation is very long 
while dynamical effects involving the SO coupling remain 
in a fixed relation to those induced by direct phonon cou- 
pling and are negligible in comparison to them. 
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